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1. Introduction 

The localization of a non-interacting wave form in a disordered potential, predicted by 
Anderson [T], has been the topic of vigorous research in different areas. Localization has 
been observed experimentally in diverse contexts [2] including Bose-Einstein condensates 
(BEC) [Hill]. Billy et al. [3] demonstrated the localization of a cigar-shaped interacting 
87 Rb BEC released into a one-dimensional (ID) waveguide with controlled disorder 
created by a laser speckle. Roati et al. [1] observed the localization of a non- 
interacting 39 K BEC in a bichromatic quasi-periodic optical-lattice (OL) potential 
created by superposing of two standing-wave polarized laser beams with incommensurate 
wavelengths. The non-interacting BEC was created [1] by tuning the atomic scattering 
length a to zero near a Feshbach resonance [5]. The disorder in a quasi-periodic OL 
potential is deterministic, in contrast to the complete disorder in a optical speckle 
potential. The localization in such a quasi-periodic potential is a special case of 
Anderson localization in a fully disordered potential and is well described by the Aubry- 
Andre model [6] . Such a localization is often termed Aubry- Andre localization. However, 
these two mechanisms of localization are distinct. While Anderson localization of wave 
functions with exponential tails is a pure quantum effect, the Aubry- Andre localization 
may occur in a classical phase space [7]. The localization of a BEC in a disordered 
potential has been the subject matter of various theoretical [El El QUI [HJ EE] and 
experimental [2j [3j studies. 

In the presence of strong disorder one has strong Anderson localization [3111], where 
the localized state could be quite similar to a localized state of Gaussian shape in an 
infinite potential or a potential of very high barriers. Then the quantum state cannot 
escape the strong barriers of the disordered potential. However, the more interesting 
case of localization is in the presence of a weak disorder when the system is localized 
due to the quasi-periodic (disordered) nature of the potential [31 H] and not due to 
the strength of the lattice. The localization takes place due to cancellation of waves 
coming after multiple scattering from many barriers of the random potential. When 
this happens the localized state acquires a pronounced exponential tail. 

The usual dilute BEC with negligible dipole moment interacting via short-range 
interaction is described by the mean-field Gross-Pitaevskii (GP) equation. More 
recently, it has been possible [13] to obtain a BEC of 52 Cr atoms with large dipole 
moment leading to a long-range dipole-dipole (dipolar) interaction superposed on the 
usual short-range interaction, which can be varied using a Feshbach resonance [5]. This 
allows to study the dipolar BEC (DBEC) of 52 Cr atoms with variable short-range 
interaction [13J. Because of the anisotropic dipolar interaction, the DBEC possesses 
many distinct features [131 [111 UHl [161 EE], which are under active investigation by 
different research groups [HI [19]. For example, the stability of a DBEC depends not 
only on the value of the scattering length, but also strongly on the geometry of the 
trapping potential [HI HU H5] . A pancake-shaped trapping potential give a repulsive 
mean field of dipolar interaction and thus the dipolar condensate is more stable. In 
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contrast, a cigar-shaped trapping potential yields an attractive mean-field of dipolar 
interaction and hence leading to a dipolar collapse [131 EOj [2H [22j [23]. Peculiarities 
in the collective low energy shape oscillations of DBEC have been studied by Yi and 
You [21]. Analogue of roton-maxon instability and the appearance of roton minimum 
in Bogoliubov spectrum [2^1 [25] and angular collapse of DBEC [26] have been studied. 
Numerical studies on the stability shows certain unusual structure of the Hartree ground 
state of dipolar condensate in an anisotropic trap [151 [16] - 

After the experiment [I] on the localization of a ID BEC in a quasi-periodic trap, a 
natural extension of this phenomenon would be to achieve localization in two and three 
dimensions (3D) [271 [T2] . both for a BEC and a DBEC. The theoretical description of 
a dilute weakly-interacting DBEC can be formulated by including a dipolar interaction 
term in the GP equation maintaining the formal simplicity, nevertheless, increasing 
vastly the numerical complexity. Using the numerical and variational solutions of the 
GP equation, we study the localization of a BEC and a DBEC in 3D, in the presence 
of bichromatic OL potentials along orthogonal directions. Although the localized states 
have an exponential tail in a weak quasi-periodic potential, the central part of such 
localized states, responsible for the major contribution to the total density, may have 
a Gaussian distribution. In the present paper we shall consider such localized states 
with an exponential tail and a Gaussian central part, which can also be studied by 
the variational approximation. The variational approximation provides an analytical 
understanding of the localization and also yields interesting result when the numerical 
procedure is difficult to implement. Such a variational Gaussian approximation has 
successfully been applied to the localization of a BEC without dipolar interaction in one 
[10] . two and three dimensions [12] in a bichromatic OL potential as well as a speckle 
potential [TT] in one dimension. 

Due to the angle dependence of the long-range dipolar interaction, the localization 
of a DBEC is more interesting than a BEC with only a short-range interaction. The 
dipolar interaction is weak in the spherically-symmetric shape, is attractive in the cigar 
shape (aligned dipoles arranged linearly attract each other) and repulsive in the pancake 
shape (aligned dipoles arranged side-by-side repeal each other) [19]. Because of this, the 
trap configuration plays an essential role in the localization of a DBEC. The controllable 
short-range interaction together with the exotic dipolar interaction makes the DBEC 
an attractive system for experimental localization in a bichromatic OL potential and a 
challenging system for theoretical investigation allowing to study the interplay between 
the dipolar interaction and the short-range interaction in dipolar atoms. 

The cigar- and pancake-shaped localized DBEC are obtained by considering 
bichromatic OL potentials of different strengths along the orthogonal directions. A 
localized DBEC, without short-range interaction, can be achieved in a bichromatic 
OL trap for a small number of atoms by tuning [13] the atomic scattering length to 
zero near a Feshbach resonance [5]. For a cigar-shaped DBEC, without short-range 
interaction, as the number of atoms increases it becomes highly attractive and suffers 
from collapse instability. In the presence of a repulsive short-range interaction, the effect 
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of the attractive dipolar interaction in the cigar shape can be compensated leading to 
a localized DBEC for a small number of atoms; derealization may take place due to 
excess of repulsion for a large number of atoms. For a pancake-shaped DBEC, the 
dipolar interaction is repulsive and a localized DBEC is obtained for a small number of 
atoms. For a large number of atoms excess of repulsion should lead to derealization. 
From a variational analysis of the localization of a DBEC, a phase diagram illustrating 
its stability for different short-range interaction, number of atoms, and the geometry of 
the bichromatic OL is given. To obtain localization, the short-range interaction should 
be small. As in ID [9j [TOj [28] , a large repulsive short-range interaction destroys the 
localization and the localized state escapes to infinity in all cases. A large attractive 
short-range interaction leads to collapse instability and destroys the localization. 

In Sec. II we present a brief account of the modified GP equation with the 
bichromatic OL potential in the presence of a dipolar interaction together with a 
Gaussian variational analysis. In Sec. Ill we present numerical and variational studies of 
localization of a dipolar dipolar interaction in the presence and absence of a short-range 
interaction. In Sec. IV we present a brief summary of the present investigation. 



2. Analytical Formulation 



We shall study the localization of the DBEC of N atoms, each of mass m, using the 
following mean-field GP equation with bichromatic OL potential V(r): [T3] 



with 



and 
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V(r) + 4vraA|0(r)| 2 + J U dd (r - r')|0(r')| 2 dr' 



^(r) = 3Aa dd (l-3cos 2 #)/r 3 , 



V{r) = Ai[v P {wc? (hx) + sin 2 (hy)} + v z sin 2 (k\Z)\, 
i=i 



(1) 



(2) 



(3) 



where <f>(r) is the DBEC wave function of normalization / (f>(r) 2 dr — 1, \i is the chemical 
potential, a the atomic scattering length, 9 is the angle between r and the direction 
of polarization, here taken along the z axis, A\ = kfsi/2, A;'s are the wave lengths, 
ki = 2n/\i are the wave numbers, and si are the strengths of the OL potentials. The 
parameters v p and v z control the relative strengths of the OL's in different directions and 
can be varied to achieve the pancake- {y z > 1, v p = 1) and cigar-shaped [y p > l,u z = 1) 
DBEC. The constant a d d = ^ojj 2 m/(12irh 2 ) in the dipolar interaction Udd is a length 
characterizing the strength of dipolar interaction and its experimental value for 52 Cr 
atoms is 15ao [13], with ao the Bohr radius, Jx the (magnetic) dipole moment of a single 
atom, and /io the permeability of free space. 

The GP equation ([!]) can be solved variationally by minimizing the energy 
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E[4>}= J -|V0| 2 + ]TV(r)0 2 + 27raA^ 4 + |- J U dd (v - v')W)\ 2 dv 

3=1 



dr 



with the following Gaussian ansatz for 0(r) [29] : 
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where w p is the width in the radial direction p and w z is the width in the axial direction 
z. We have assumed circular symmetry in the x — y plane. Equations fll]) and (jSJ) lead 
to [29] 
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where £ P: i = exp(— w 2 kf), £ z ,i = exp(— u> 2 /cf ). For a stationary state the energy mE/h 
should have a minimum as a function of w p and w z : dE/dw p = dE/dw z = 0, together 
with the condition 
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The minima conditions dE/dw p = dE/dw z = become in explicit notation 
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The solution of © - 
parameters, localization is possible if there is a minimum of the energy. Hence from 
and — (fT2"|) we can determine if a BEC or DBEC would be localized or not. 



(12) 



12"]) determine the widths and For a certain set of 



3. Numerical Study 



We perform a full 3D numerical simulation in Cartesian x, y, z variables using imaginary- 
and real-time propagation with Crank-Nicolson discretization [30] employing small space 




Figure 1. (Color online) (a) Normalized numerical (Num) and variational (Var) 
densities (j) 2 (x) of (fT3|) versus x for the ID localized state with bichromatic OL wave 
lengths Ax = 5, A 2 = 0.862A;l, and strengths S\ = s 2 = 4. (b) The same densities 
together with the exponential fit <j) 2 {x) = cxp[— 2abs(a;)/Zi oc ]/(0.4865v / 7r), hoc = 0.75 
on a log scale. 

(~ 0.025) and time (~ 0.0001) steps necessary for obtaining converged results. For this 
purpose we use the FORTRAN programs provided in [31] after transforming (CQ) into 
a time-dependent form by replacing m/i/h 2 by id/dt, where t is time. The dipolar 
interaction term is evaluated by the usual fast Fourier transformation technique (14] . 
The imaginary- and real-time propagation lead essentially to the same localized states. 
This not only assures the correctness of the calculational scheme, but also guarantees 
that the localized states are stationary and not a result of dynamical self-trapping 
(32] . (The imaginary-time propagation can only find the stationary localized states and 
cannot obtain the dynamical self-trapped states.) The stability of the localized state was 
tested by real-time propagation allowing small perturbations of potential or interaction 
parameters. (In the absence of the dipolar interaction, the localized states have been 
demonstrated explicitly to be stable in one \TU\ and two p2] dimensions.) The accuracy 
of the numerical simulation was tested by varying the size of space and time steps and 
the total number of space and time steps. 

To compare the numerical results with variational analysis, we only consider 
localized states mostly occupying a single site of the bichromatic OL. Throughout this 
study the strength parameters of the two components of the bichromatic OL are taken as 
Si = s 2 = 4 and the corresponding wave lengths are taken as Ai = 5 fim, X 2 = 0.862Ai. 
The relative strength of the radial and axial bichromatic optical lattice are varied by 
adjusting the parameters v p and v z in Eq. ([3]). To demonstrate that with these sets 
of parameters we are in the limit of localization with an exponential tail in a weak 
potential, we solve the ID linear Schrodinger equation in dimensionless variables [31] 

- + E ^ sin 2 C^x) <f>(x) = £<f>(x), (13) 

with the above sets of parameter, where the prime denotes x-derivative and £ denotes 
energy. In the limit of zero nonlinearity (a = 0) and zero dipolar interaction (a^d = 0), 
([I]) decouples into three equations like ( TT3"|) . In figure [TJ (a) and (b) we plot the density 
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Figure 2. (Color online) Variational widths w p and w z of the DBEC versus 
bichromatic OL relative strength parameters v p and v z [viz. Eq. (j^J] for 170 52 Cr 
atoms for scattering length a = and 20ao. 

<ft 2 (x) of ([TBI versus x in linear and log scales together with its Gaussian variational 
counterpart <p 2 (x) = exp(—x 2 /w 2 )/(7i 1 ^ 2 w), with w = 0.4865 the variational width, 
and the exponential fit <p 2 (x) = exp [-abs(a;)//i oc ]/(0.48657r 1 / 2 ), with l loc = 0.75 the 
localization length [3] providing a measure of the exponential tail. For strong localization 
with strong bichromatic potential l\ oc — > and the localized state has a pure Gaussian 
tail. However, when l\ oc > x rms with x rms the root mean square size, the exponential tail 
is pronounced and the limit of localization in a weak quasi-periodic potential is attained 
[3111]. In the present example, x rms = 0.53 and Zi oc = 0.75, and hence we are in the limit 
of localization in a weak potential. This is clear from figure [H (b), where the central 
part of the localized state is fitted to the variational Gaussian solution, whereas the 
large-x parts are fitted to an exponential function with a large localization length over 
about nine orders of magnitude. In figured] (b), one can identify several minor peaks in 
successive wells of the bichromatic OL potential. The state of figured] is localized by the 
weak bichromatic lattice and if we substantially reduce the strength of the potential, no 
localized state will emerge. In the experiments of Refs. [3] 0] and in some other studies 
on localization [9] in a weak potential, localized states with a pronounced undulating tail 
over many wells of the localization potential were considered. However, a weak limit of 
localization with an exponential tail can also be achieved in the absence of a pronounced 
undulating tail, as we have shown here. In this paper we shall be considering localization 
in the presence of a pronounced Gaussian peak and a weak undulating tail as can be 
seen from figured] (b). The existence of a pronounced Gaussian peak will be turned to 
a good advantage in predicting accurate analytical variational results for the localized 
state and for an analytical understanding of the localization. We shall also consider 
the localization in the presence of short-range interaction and dipolar interaction. The 
inclusion of a repulsive short-range interaction will in general destroy localization by 
increasing the localization length and thus creating a more pronounced exponential tail 




Figure 3. (Color online) Stability region in the a/ao versus v p and v z (relative 
strengths of the bichromatic OL in radial and axial directions) phase plots for (a) 170 
and (b) 500 atoms for add = 15ao ( 52 Cr atoms) and add = variational - var (lines), 
numerical - num (points). Localization is possible between the upper and lower lines 
of the same data set. 



[9] . The inclusion of the dipolar interaction will have an effect on localization, which we 
shall study here, and such inclusion should not destroy the exponential tail of localization 
as illustrated in figured] (b). 

A variational analysis is useful for a qualitative understanding of the problem and 
we present the same before considering a numerical solution of (CQ). In figure [2] we plot 
the variational widths w p and w z for a DBEC of 170 52 Cr atoms in the cigar and pancake 
shapes. This figure shows the evolution of the widths from the pancake to cigar shapes 
while the radial width w p reduces and the axial width w z increases as expected. For 
a = 0, the axial width w z does not increase with the increase of u p , as in the cigar shape 
the dipolar interaction becomes attractive and does not permit the increase of w z . 

Next we study, using flSJ), the set of values of the parameters for which the energy 
can have a minimum and allows a stable localized DBEC. The region of stability for 
dipole strength add = and = 15ao is shown in figure [3] for (a) 170 and (b) 500 52 Cr 
atoms as a phase plot of a/ao versus the relative strengths of the bichromatic OL v p 
and u g in radial and axial directions. The stability of a DBEC in a harmonic trap has 
also been studied [H]. For a fixed add and N the stability region lies between the two 
corresponding lines in figure [3j Above the upper line the system becomes too repulsive 
(positive Na/ao) to be confined by the weak bichromatic OL. Below the lower line the 
system becomes too attractive (negative Na/ao) an d suffers from collapse instability. 
The localization of a BEC (without dipolar interaction) is controlled by the non-linearity 
AnaN alone of ([T|) and the effect of the dipolar interaction on the stability of the DBEC 
is clearly exhibited in figures [3j In these figures we also plot numerical results for collapse 
instability with negative (attractive) scattering length (the lower limit of stability in this 
figure), which agree well with the variational results. Starting from a stable localized 
state, the stability lines (points) are obtained by slowly changing the scattering length 




Figure 4. (Color online) Numerical (num) and variational (var) energy E versus 
relative strengths of bichromatic OL v p and v z for (a) 170 and (b) 50 52 Cr atoms. 

a at a fixed trap symmetry (by fixing the parameters v p and v z ) until no localized 
state can be obtained (by numerical or variational means). In the pancake side the 
dipolar interaction is repulsive and the localization is destroyed for a smaller value of 
the repulsive scattering length compared to the case where dipolar interaction is absent, 
as can be seen in figure EJ The opposite happens in the cigar side where the dipolar 
interaction is attractive. In the cigar side the dipolar interaction is attractive and the 
localization is destroyed for a larger value of the repulsive scattering length compared 
to the case where dipolar interaction is absent, as can be seen in figure |3] (b) for 500 
atoms. This effect is smaller in figure [3] (a) for 170 atoms. 

From (ED) we find that f(n) is positive for cigar shape leading to an attractive 
contribution to energy fl6]), whereas it is negative in the pancake shape leading to a 
repulsive term in energy. For N = 170 there is a symmetric state with w z = w p and 
f(n) = near v p — « 1.3 where the dipolar interaction does not contribute and 
where the DBEC is most stable at the maxima in figure |3] (b). At this point the DBEC 
acts like a "normal BEC" and the Na/ao value at the maximum is independent of N. 
In the pancake shape, the localization of the DBEC can be easily destroyed due to a 
large repulsive dipolar interaction in the weak bichromatic OL. In the cigar shape, the 
large attractive dipolar interaction may lead to the collapse of the localized DBEC more 
easily than in the absence of dipolar interaction. These aspects are clearly illustrated 
in figure |3] (b). For iV = 500 the interactions are stronger than for N = 170 and the 
domain of allowed localization in terms of number of atoms has reduced. 

In figures |4] (a) and (b) we compare the numerical and variational energies of the 
52 Cr DBEC of 170 atoms with a = and 20ao and of 500 atoms with a = and 10ao, 
respectively. A smaller value of a is required for the stability of 500 atoms [see figure 
|3] (b)]. The agreement between variational and numerical results is good in general 
except for v z — > 100, v p = 1. In this limit the localized DBEC occupies two sites of the 
bichromatic OL and does not have a Gaussian shape. This justifies a small disagreement 
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Figure 5. (Color online) 2D Contour plot of density 4> 2 (x, 0, z) in the y = plane for 
the cigar-shaped DBEC with v p = 10, v z = 1 of N = 500 52 Cr atoms with a — 0: (a) 
Numerical (b) variational. 2D Contour plot of density <f> (x,y, 0) in the z = plane 
for the same DBEC: (c) Numerical (d) variational. 

between the numerical and variational results in this limit. The variational result for 
a = 10a in figure @] (b) only exists up to v z pa 35 [see figure |3] (b)]. 

Next we study how the DBEC solely under the effect of dipolar interaction (a = 0) 
changes its shape as we move from the cigar to pancake shaped configuration. In 
figures [5] (a) and (b) we show the 2D contour plot from the numerical and variational 
analysis, respectively, for density 4> 2 (x, 0, z) in the y — plane for a cigar-shaped DBEC 
with Up = 10, v z = 1 for 500 52 Cr atoms for a = 0. The 2D contour plot from the 
numerical and variational analysis, respectively, for density <p 2 (x,y,0) in the z = 
plane for the same DBEC is shown in figures [5] (c) and (d), respectively. This localized 
numerical DBEC state is of small size, compared to its variational counterpart, due 
to the attractive dipolar interaction in the cigar-shaped DBEC. The attraction due to 
dipolar interaction has caused the DBEC to contract from an average Gaussian shape. 
The corresponding numerical and variational energies of the localized states of figures 
[5] are 15.16 and 15.53, respectively. The localized DBEC states in figures [5] (and also in 
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Figure 6. (Color online) Same as in figure [5] for the DBEC in the pancake trap with 
v p = l,v z = 10. 

figures [6]) occupy practically a single bichromatic OL site at the origin along x, y and z: 
—1.25 /iin < x,y, z < 1.25 /xm. However, these localized states have exponential tails. 

The effect of the dipolar interaction is small in the symmetric case with v z — v p — 1 
and N = 500. Nevertheless, due to the complicated dipolar interaction, the density 
<j) 2 (x, 0, z) does not have a pure circular shape in this case (not explicitly shown). This 
density distribution will be circular for a BEC without dipolar interaction. (A nearly 
circular shape for this density is obtained for a DBEC with v p = 1 and v z = 1.6.) The 
numerical and variational energies for this state are 6.97 and 7.06, respectively. 

We next consider the localization for v z = 10, v p = 1 and a = 0. In this pancake 
shape, the dipolar interaction is repulsive and for 500 52 Cr atoms a localized DBEC is 
obtained. In figures [6] (a) and (b) we show the numerical and variational 2D contour 
plots of density 4> 2 (x, 0, z) in the y = plane, respectively, for this DBEC. The numerical 
and variational 2D contour plots for density 2 (x, y,0) in the z = plane for the 
same DBEC is shown in figures [6] (c) and (d), respectively. In this case, the dipolar 
interaction is repulsive in nature and hence the size of the DBEC is larger than the 
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Gaussian variational shape (see figures \5\ where the opposite happens for a DBEC with 
attractive dipolar interaction.). The numerical and variational energies for this state 
are 14.16 and 14.26, respectively. The change from a cigar shape to a pancake shape 
of the localized DBEC as we move from v p = 10 and v z — 1 to v p — 1 and u z = 10 is 
obvious from the density distribution in figures |5] (a) and [6] (a). Although the variational 
and numerical energies (studied in figures |5] and [6]) are quite close to each other, the 
numerically obtained matter density should have some peculiarities not obtainable from 
the variational calculation due to the anisotropic dipolar interaction. (The variational 
calculation is based on an axially-symmetric Gaussian distribution.) 

4. Summary and Conclusion 

We investigated the localization of a 52 Cr DBEC with in a weak bichromatic OL trap in 
the presence and absence of short-range interaction using the numerical and variational 
solution of the 3D GP equation ([T|). Of the two solutions, the numerical solution is 
the most precise one and should be used in case of disagreement with the variational 
solution. Although the density of the central part of the localized states has a near 
Gaussian distribution, the density distribution also has a long exponential tail [3j 0]. 
The Gaussian distribution near the center permits a variational analysis of localization, 
which is used for an analytical understanding of the problem. A DBEC of a small number 
of atoms with a weak short-range interaction could be localized by a relatively weak 
bichromatic OL trap. From the variational solution we obtain a phase diagram [figure 
|3] (b)] illustrating the effect of the dipolar interaction on the localization as a function 
of the strengths of the trap u z and v p in axial and radial directions, respectively. We 
find that for 52 Cr atoms, the dipolar interaction has a moderate effect on localization. 
(Larger effect will certainly appear for dipolar molecules where the dipole moment could 
be larger by an order of magnitude compared to the dipole moment of 52 Cr atoms.) The 
numerical and variational energies of the DBEC, as well as the corresponding densities 
are in reasonable agreement with each other. In the absence of a short-range interaction, 
the localized DBEC can accommodate the largest number of 52 Cr atoms (~ 1000) in the 
spherical configuration and this number reduces both for cigar and pancake shapes due 
to the attractive and repulsive dipolar interaction, respectively. The attractive dipolar 
interaction leads to collapse and the repulsive dipolar interaction leads to leakage to 
infinity. We hope that this study will motivate experiments on the localization of a 52 Cr 
DBEC in a bichromatic OL trap. The estimate of the number of localized 52 Cr atoms, 
their radial and axial sizes and shapes etc. as predicted in the present study can be 
verified in the experiment. 
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